The struggle for space: Viral extinction through competition for cells 
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The design of protocols to suppress the propagation of viral infections is an enduring enterprise, especially 
hindered by limited knowledge of the mechanisms through which extinction of infection propagation comes 
about. We here report on a mechanism causing extinction of a propagating infection due to intraspecific compe- 
tition to infect susceptible hosts. Beneficial mutations allow the pathogen to increase the production of progeny, 
while the host cell is allowed to develop defenses against infection. When the number of susceptible cells is 
unlimited, a feedback runaway co-evolution between host resistance and progeny production occurs. However, 
physical space limits the advantage that the virus can obtain from increasing offspring numbers, thus infection 
clearance may result from an increase in host defenses beyond a finite threshold. Our results might be relevant 
to better understand propagation of viral infections in tissues with mobility constraints, and the implications that 
environments with different geometrical properties might have in devising control strategies. 

PACS numbers: 87.23.-n, 87.10.-e, 87.18.-h 



One of the common strategies used by viruses to ensure 
fast adaptation is the steady production of mutants and of 
enormous amounts of progeny, two characteristics that chal- 
lenge the design of therapies able to induce viral extinction. 
Increased mutagenesis has been successfully used in vitro to 
eradicate infectivity, though there is no general agreement on 
the pathways through which extinction supervenes |[T}21. In 
mutational meltdown, the population becomes extinct because 
the average growth rate falls below one, through a process 
different from the error threshold predicted by classical qua- 
sispecies models |5 1 or from stochastic extinction of infectiv- 
ity mi?]- The progress of an infection is further conditioned 
by the geometry of the space where it occurs lIS HTOl . Experi- 
ments with viruses infecting bacteria have highlighted differ- 
ences in the internal structure of populations evolving in liquid 
medium or on bacterial monolayers ifTTI . These results might 
give clues to treat viral infections of two-dimensional tissues: 
It has been shown that space induces a remarkably strict clus- 
tering of the propagation in crops lfT2l . in leaves lfT3l . and in 
lawns of cells [T?]. Advances in therapy design ask for an 
improved understanding of the mechanisms causing viral ex- 
tinction, before they can be applied in vivo. 

Models of viral evolution necessarily make simplifying as- 
sumptions, and real virus behavior often deviates substantially 
from their predictions lITSl . Current quasispecies models im- 
plement high mutation rates that yield heterogeneous popula- 
tions — a property of most viruses infecting plants lfT6l and 
of many viral populations on cellular monolayers fYf\. One 
common approximation is to assume that all new mutations 
have a deleterious effect on fitness, thus neglecting beneficial 
and neutral mutations. However, adaptation occurs frequently, 
and even at low population numbers. Classical quasispecies 
models assumed the existence of a unique master sequence 
of high fitness, while we cuiTently know that a huge amount 
of different sequences yield phenotypes that perform equally 
well. When beneficial and neutral mutations are considered, 
extinction might supervene through mechanisms other than 



error thresholds [4]. The predictions yielded by models of 
viral evolution depend on whether physical space is explicitly 
considered or not |18|. In quasispecies models, cell-to-cell 
transmission or local transport of viruses unavoidably results 
in a lower value of the critical error threshold liT9ll2()l . Clus- 
terization of viral types is responsible for the latter effect, as 
well as for enhancing the coexistence of types in the quasis- 
pecies iioiisn . 

In this Letter we introduce a model representing the prop- 
agation of infections on two dimensional arrays of cells. We 
incorporate two realistic features of infection propagation of- 
ten disregarded in quasispecies models: a non-negligible frac- 
tion of beneficial mutations and a host resistance to infection. 
We show that with spare susceptible cells at each generation, 
the viiTis can escape host resistance by increasing its progeny. 
However, when viruses compete for infecting cells, increasing 
the progeny does not pay and infection clearance occurs at a 
finite value of the host resistance. In this case the model boils 
down to a multicomponent generalization of the Domany- 
Kinzel (DK) probabilistic cellular automaton ll22l hence clas- 
sifying viral extinction within the directed percolation (DP) 
universality class, alongside with propagation of immiscible 
fluids in porous media and other similar problems. 

In the model, a viral population is described as an ensem- 
ble of (pheno)types, each characterized by a replicative abil- 
ity r = 0, i? standing for the number of offspring able to 
infect cells produced by each type under replication. Delete- 
rious or beneficial mutations decrease or increase (with prob- 
abilities p and q respectively) offspring's replicative ability in 
one unit. Susceptible cells have a resistance to infection tt, 
i.e., a viral particle can infect a cell with probability 1 — tt. 
For simplicity, we consider that the dynamics proceeds in dis- 
crete generations. 

Assume first that the availability of susceptible cells is large 
enough for every viral particle with r > to actually meet 
a cell to infect, no matter its replicative ability r. This is 
a mean-field dynamical description, where infection occurs 
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with a probability 1 — tt. If n,. (t) denotes the number of viral 
particles of type r at generation t, for 1 < r < R 

nr{t + 1) = (1 - tt) [r(l -p- q)nr{t) 

+ (r + l)prv+i(0 + (r - l)(7nr_i(i)] , 
nR{t + 1) = (1 - tt) - p)nR{t) + {R- l)qnR^i{t)] . 

(1) 

Besides, there is a class r = which has lost its replicative 
ability and is maintained by class r — \ through deleterious 
mutations, i.e. ?T.o(t + l) = (1 — 7r)pni(i). As initial condition 
we take nr{Q) = Sr.R- 

An exact analytic solution exists for the case q = which 
represents the worst situation for the survival of the virus: if 
extinction does not occur for g = 0, it will not occur for g > 
either. In the limit t — ^ oo, the population of every class grows 
at a fixed rate, 71^(^ + 1) = \nr{t). System ([TJ then becomes 
an eigenvalue problem, A being its largest eigenvalue. For 
q — Q, the asymptotic growth rate is A = (1 — 7r)i?(l ~p) and 
the fraction of virus in class r is 



= lim 



nr{t) 



t->-oo 



R 



pR-^l^py. (2) 



Extinction occurs when A < 1, i.e., if the maximum progeny 
production R < R^ = (l-7r)~i(l-p)~^ Although the 
larger the resistance to infection of the cell tt the more de- 
manding becomes this condition, R can always be increased 
beyond Rc to avoid extinction. If the number of suscepti- 
ble cells is unlimited and increasing R always improves the 
chances of survival of the pathogen, a runaway co-evolution 
between the virus and the host is to be expected. 

The case q ~ is not representative of the situation when 
q > though. The stationary distribution (|2]i has a well de- 
fined R oo limit, whereas the permanent improvement al- 
lowed when g > forbids a steady state in the same limit. 
As a matter of fact, it can be shown that in that case the viral 
population grows superexponentially |23|. 

Consider now a restricted number of susceptible cells avail- 
able. This occurs when cells have a spatial distribution and/or 
when the mobility of viruses is limited. To assess the con- 
sequences of this constraint, we study the dynamics proceed- 
ing on a triangular lattice with periodic boundary conditions. 
Even if infection starts off from a single infected cell, the 
propagating front becomes asymptotically flat LlOl . Hence, 
we consider that cells form arrays in one spatial dimension, 
and infection advances one row of cells per generation. Let 
ri (t) be the replicative ability of the individual occupying site 
j = 1, . . . , L at generation t. 

To account for the limited mobility of the viral particles, 
they are constrained to infect only neighboring cells. Hence, 
at generation t + 1 site i will be occupied by one of the off- 
spring of the individuals at sites i or i + 1 at the parental 
generation t. The probability that ri{t + 1) — r, given that 



ri{t) = ri and n+i{t) = r2, is 

P{r\r,,r2) o + (l - 7r"^+"^) ^ "^^^""^ , 

' ^ ^ —I— r-r-. 



»'i + r2 
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where Pr,r-i = P for f < R, Pr.r+i = q and p^.r = 1 ^ P ^ 9 
for r < R,pR^R+i = 0, andp^.^, = 1—p. We take ri(0) = R 
for all i — 1, . . . ,L when the infection begins. According 
to ([3|, the probability that the cell resists the viral attack is 



given by vr 



ri+r2 



i.e. the probability that infection fails after 



ri + r2 independent trials. With the complementary prob- 



ability 1 



infection occurs, and one individual ran- 



domly drawn among the ri + r2 possible occupies the site. Its 
replicative ability may change due to mutations as specified. 
The limit case n ^ 1 describes an immune host, whereas in 
the limit tt — > the host plays no role in the progress of infec- 
tion. Two relevant quantities that characterize the dynamical 
state of the system are the spatial density p{t) of virus with 
r > 1 at time t and the average replicative ability of the pop- 
ulation f{t) = {p{t)L)^^ A systematic study of 
these two quantities in the limit i — > oo as a function of p and 
TT portrays the asymptotic behavior of the model. 

In the absence of beneficial mutations (q — 0), the class 
with the highest r present disappears when p or tt increase be- 
yond a threshold value. These transitions are somehow anal- 
ogous to the error threshold described for simple models of 
quasispecies |5|, where the master sequence disappears and 
the population turns into a cloud of low-fitness replicating mu- 
tants. For constant resistance to infection tt, there is a hierar- 
chy of domains in the phase space where high— r classes are 
sequentially lost as p increases (see Fig.[T]i. For increasing tt, 
the largest r— value present at global extinction increases. 

There are two features of Fig. [T] for q = which are 
amenable to further analysis. The first one is the last error 
threshold, when class r = 1 disappears. There is a mapping 
of this transition to the DK cellular automaton. When only 
classes r — 1 and r = are present, sites are either replica- 
tive (occupied) or non-replicative (empty), yielding probabil- 
ities ^ P(1|0,0) = 0, pi = P(1|1,0) = P(1|0,1) = 
(l-7r)(l-p), andp2 = 1) = (1 - 7r2)(l -p), which 

define a DK automaton. Curve C in Fig. [TJcoincides with the 
DK transition cast in the variables (p, tt) describing our model. 
This curve meets the tt = axis at p^ = 0.294510(6), a value 
derived from the known site percolation critical point |24 |. 

The second feature obtained from a mapping with the DK 
model is the extinction threshold of the highest class r = i? at 
p = 0. Since the initial state only contains individuals of this 
type, no other types will ever appear because of the lack of 
deleterious mutations. Hence, sites can be either occupied or 
empty, just as in the DK model. Transition probabilities corre- 
spond to the bond percolation case, from which it follows that 
the thi-eshold occurs at tTc = a^^'^, with a = 0.3552998(9). 
Only in the limit i? — > oo would tTc 1, taking a lower value 
for any other combination of parameters. 

For p and g > 0, no matter how small, all replicative 
types are simultaneously present in the infective phase (see 
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FIG. 1. (Color online) Phase space of the model with R — 10 and 
two choices for q: q — 0, where the highest— r classes are sequen- 
tially lost as p increases through transitions analogous to the error 
threshold of Eigen's quasispecies model (white lines), and q = 0.01, 
where only global extinction of the population is possible. Colors 
code for the average replicative ability f when g = 0.01. In the 
white domain, infection is cleared in finite time. The top extinction 
occurring for q = can be obtained from that of the DK cellular 
automaton (curve C). The inset shows f for p — 0.11 and two val- 
ues of q. Data have been obtained with a system of size L — 10* 
and averaging over 5 x 10'^ time steps after a transient of 2 x lO** 
generations. 

Fig.[T]i. For a fixed value of vr, lower r classes get more pop- 
ulated upon increasing p, hence the average replicative ability 
decreases. The behavior as tt increases at constant p is of a dif- 
ferent nature. As tt increases, infected cells become sparser. 
In the absence of neighbors, the probability that a cell at gen- 
eration t + 1 is infected by offspring of an individual of class 
r at generation t is proportional to (1 — tt*"). This quantity pe- 
nalizes more severely low— r classes as tt increases, such that 
high— r classes receive a relative advantage and f increases 
as TTc is approached. At the critical value tTc the quasispecies 
eventually collapses because propagation is prevented by host 
defenses (see the inset of Fig.[T]). 

There is another relevant difference between the two mod- 
els. As we have seen, when cells are in excess any host re- 
sistance TT can be circumvented by increasing R. In contrast, 
when viral classes have to compete for susceptible cells, there 
is a critical value of tt above which this is no longer true. The 
reason (illustrated by Fig. |2]i is that the average replicative 
ability f quickly saturates to a constant value upon increas- 
ing R. In other words, when viruses struggle for infecting a 
limited number of cells, increasing R does not provide any 
additional advantage. In order to show that this effect results 
from the limitation in the number of cells available, we have 
modified the model to consider a case where four cells are 
available for each infection event. This amounts to increasing 
the competition between types, such that low— r classes are 
at a disadvantage with respect to the previous situation and f 
increases ifTOl (see Fig.[2]i. As a side effect, the critical value 
of the cell resistance, tTc, at which extinction occurs also satu- 
rates to a value below 1 as i? oo. A runaway co-evolution 




FIG. 2. When the number of susceptible cells available is limited, 
the average replicative ability f (open symbols, left j/— axis) and the 
critical value of the resistance to infection ttc (solid symbols, right 
y— axis) saturate to finite values as 7? increases. In these curves, 
p = 0.2, q = 0.01, and averages over 50 independent realizations 
for each value of ttc have been performed. Mean-square deviations 
are smaller than symbol size. 



is therefore excluded in this case: if the host cell is able to in- 
crease its defenses beyond tTj, the infection will not progress. 

When the number of susceptible cells is limited, the tran- 
sitions to extinction belong to the universality class of DP for 
any q > 0. For q = we have rigorously proven this for 
the extinction of the r — 1 class as well as when p — > 0, 
by mapping the model to a DK cellular automaton, known to 
belong to the DP class ll24]| . Besides, our model is defined 
on a set of local rules, has no special symmetries or disor- 
der, and exhibits a continuous transition to a unique absorb- 
ing state (total extinction). Therefore, according to the DP 
conjecture it should belong to the DP universality class. An 
empirical proof is provided in Fig. [3] where p{t) is appropri- 
ately rescaled according to the exponents of DP. The transition 
point pc is numerically determined as the value of p yielding 
the best collapse. Finally, the study of the variation of p{t) 
at Pc with the system size L is characterized by the dynamic 
exponent z — 1.580745(10). This analysis has been repeated 
for the transition where r — 2 disappears at q = 0, with anal- 
ogous results. 

Formally, the probabilistic model presented here defines a 
multicomponent generalization of the DK cellular automaton. 
Our model only focuses in a limited region of its parameters in 
which an active-inactive transition (q > 0) or several (q = 0) 
occur with the criticality of directed percolation. Exploring 
in depth the whole phase space of this model in search for 
deviations of this behavior is an interesting open question left 
out by the present study. 

Our model could be modified to better mimic other biolog- 
ical situations. In many plants, different individuals display 
variable degrees of resistance to viral infection |25|, a case 
that could be implemented by introducing a host-dependent 
TT. This would introduce a form of quenched spatial disorder 
that may lead to universality classes other than DP |24|. Anal- 
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FIG. 3. Data collapse of curves p{t) according to the critical ex- 
ponents known to characterize (1 + 1) directed percolation, which are 
5 = 0.159464(6), = 1.733847(6), and z = 1.580745(10) |24|. 
Parameters in these simulations are _R = 2, q = 0.01, tt — 0.3. The 
system size is L = 10'* for the main plot, where the legend shows 
the values of p for each curve. For this choice the transition occurs at 
Pc — 0.16600(3). Legend in the inset gives the values of the system 
size L. Data have been averaged over 50 independent realizations 
and up to 10* generations for each curve. 

ogous changes would be effected by temporally quenched dis- 
order, a situation holding, e.g., when susceptibility to infection 
depends on the age of the host. An example of the latter is 
systemic acquired resistance, a durable form of immunization 
observed in plants f26]. Spatially quenched disorder could 
change the properties of viral extinction to those of dynamic 
percolation lIZTl . Current knowledge on the properties of per- 
colation phenomena in a variety of different situations might 
inspire unprecedented strategies to stop infection propagation 
in those different environments. 

The hindrance to the viability of a virus brought about by 
a limited availability of cells has some empirical support. 
While all plant viruses infect the tissue that carries nutri- 
ents through the vascular system (phloem), where mobility is 
high and susceptible cells are abundant, some viruses (e.g. lu- 
teoviruses and some geminiviruses ll28l ) are not able to infect 
leaves, where mobility is limited to cell-to-cell transmission 
and fewer uninfected cells are available per infective cycle. 
The reverse situation is as yet unknown. 
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